Asymmetric metal-insulator transition in disordered ferromagnetic films 
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We present experimental data and a theoretical interpretation on the conductance near the metal- 
insulator transition in thin ferromagnetic Gd films of thickness b ss 2 — Wnm. A large phase 
relaxation rate caused by scattering of quasiparticles off spin wave excitations renders the dephasing 
length L$ < b in the range of sheet resistances considered, so that the effective dimension is d — 3. 
The observed approximate fractional temperature power law of the conductivity is ascribed to the 
scaling regime near the transition. The conductivity data as a function of temperature and disorder 
strength collapse on to two scaling curves for the metallic and insulating regimes. The best fit is 
obtained for a dynamical exponent z ~ 2.5 and a correlation length critical exponent u' ~ 1 .4 on 
the metallic side and a localization length exponent v w 0.8 on the insulating side. 
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Since the scaling theory of Anderson localization has 
been proposed in 1979 [![, the metal-insulator transition 
in disordered conductors [2| has been one of the most ex- 
tensively studied cases of quantum phase transition, both 
experimentally and theoretically. In its simplest form it 
describes non-interacting electrons in a disordered poten- 
tial, where the disorder can be controlled experimentally 
in a variety of ways, e.g by systematic doping. One of 
the most dramatic predictions of the scaling theory is the 
absence of extended states, and therefore true metallic 
behavior, in systems in dimensions d < 2. This has been 
verified in numerous experiments Q. The other predic- 
tion is the existence of a critical point in d > 2 where 
the conductivity in the metallic phase goes to zero con- 
tinuously with increasing disorder, in contrast to having 
a minimum metallic conductivity Electron-electron 
interactions are known to modify the behavior near a 
metal- insulator transition in a significant way Q . For ex- 
ample, indications of a metallic state in two-dimensional 
systems have been found in experiment, and a number 
of theoretical scenarios explaining such a state have been 
developed [||. 

Near the transition, the behavior is characterized by 
power laws with critical exponents. For example, the dc 
conductivity u(A), with A being a measure of disorder, 
follows a power law a ~ t s , where t = (1 — A/A c ) denotes 
the distance to the critical point at the critical disorder 
A c and s is the conductivity exponent. The dynamical 
conductivity at the critical point, on the other hand, is 
characterized by the dynamical exponent z as er(w; A c ) ~ 
a; 1 / 2 . The correlation length on the metallic side (A < 
A c ) diverges at the critical point as £ ~ t~ v and the 
localization length (A > A c ) diverges as £ ~ \t\~ V ' The 
critical exponents v and v' may be different. In d = 3 
dimensions the relation s = v 1 holds. The exponents in 
d = 3 have not been calculated in a reliable way up to 



now. 

As for any quantum phase transition, the critical ex- 
ponents can not be measured experimentally at the true 
T = critical point, but must be inferred from finite 
temperature T measurements. Therefore, the emphasis 
has been to obtain the conductivity as a function of T, 
as close to T = as possible. In spite of intense efforts 
over several decades [7|] it turned out to be rather difficult 
to access the critical regime in a reliable way. While so 
far all such experiments confirm the continuous nature 
of the transition, the values of the critical exponents re- 
main controversial. Published experimental values of s 
and z v ary from s ps 0.5 [g|, s ps 1 [9J] to s ps 1.6 and from 
z ps 2 [lj| to z m 2.94 HH- For the Anderson transi- 
tion (omitting interaction effects) , numerical studies find 
the conductivity exponent s ps 1.6 [12] while theoretical 
prediction for the dynamical exponent is z = 3 [T3 |. 

Theoretically, the scale dependent conductivity at fi- 
nite T is obtained from cr(cj) by replacing to by T. 
The critical dynamical scaling is found even away from 
the critical point, at frequencies uj > w^, where lu^ = 
r(£/0~ z - Here I and r arc the mean free path and the 
momentum relaxation time, respectively. More precisely, 
in the above scaling regime the dynamical conductivity 
obeys the scaling law, 

cr(a;;A) = r 1 G(±l ! ^ 1/z ), Q 0. (1) 

At the critical point, when £ — > oo, it follows that 
C?(±l, £ccd/ z ) ~ ^u; 1 / 2 . Using the sheet resistance i?o as 
the disorder parameter controlled in experiment, so that 
£ oc \Rq — R c \~ y , and replacing to by T, the conductivity 
should obey the scaling, 

\R - R c \- u <j(T; Ro) = G(±l, \R - R^T 1 '*) (2) 

with (Rq — R c ) ^ 0, where R c is the critical resistance. 
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Below wc will see that the exponent v can be different 
on the two sides of the transition. 

In this letter we report the study of the conductivity 
near the metal-insulator transition in a thin-film geom- 
etry where it is possible not only to increase the dis- 
order directly by changing the atomic structure, e.g. by 
varying the film deposition parameters and by annealing, 
but also indirectly by varying the film thickness. Both 
factors determine the sheet resistance, which is the sin- 
gle important parameter controlling the distance to the 
critical point. We are able to measure the disorder as 
well as the temperature dependence of the conductivity 
quite reliably and reproducibly. In particular, we are able 
to determine the location of the critical point with un- 
precedented precision. As we will explain below, at the 
relatively high temperatures and large sheet resistances 
of the experiment, the films are in an effectively three- 
dimensional regime. 

Two series of thin films of Gd (series 1 and series 2) 
were grown by r.f. magnetron sputtering through a 
shadow mask onto sapphire substrates held at a temper- 
ature of 130 K. The current and voltage leads of the de- 
posited sample overlapped with predeposited palladium 
contacts, thus allowing reliable electrical connection with 
low contact resistance for in situ measurements of the 
electrical properties. The experiments were performed in 
a specialized apparatus in which the sample can be trans- 
ferred without exposure to air from the high vacuum de- 
position chamber to an adjoining low-temperature cryo- 
stat and electrically reconnected for transport measure- 
ments. Immediately after deposition, the samples were 
transferred to the cryostat and held at a temperature of 
77 K or below. At these temperatures the samples are 
stable and do not undergo any time-dependent changes 
in resistance. If however the temperature is temporar- 
ily raised back to the deposition temperature (130 K), 
annealing marked by a slow irreversible increase in resis- 
tance occurs. 

To parameterize the amount of disorder in a given 
film[14|. we use the sheet resistance Rq = R XX (T = 5 K) 
where R xx is the longitudinal resistance. In our experi- 
ments Rq spans the range from 4 kfl (35 A thick) to 40 kO 
(< 20 A thick). Controlled thermal annealing thus allows 
us to advantageously tune a single sample through suc- 
cessive stages of increased disorder. Our series 1 samples 
comprise 5 separate depositions with two of the samples 
undergoing 12 successive anneals thus giving a total of 17 
measurements at different stages of disorder. Our series 2 
samples comprise a single sample undergoing 15 succes- 
sive anneals for a total of 16 measurements spanning the 
critical region where the metal-insulator transition oc- 
curs. The resistance was measured using four-terminal 
dc techniques. 

We now argue why the films are effectively three di- 
mensional. As discussed in Ref. jTij], the phase relax- 
ation rate in ferromagnetic films is dominated by the 



scattering off spin wave excitations, yielding r^ 1 oc T, 
in d = 2,3 dimensions. The effective dimensionality of 
the system is d = 3 if the temperature dependent correla- 
tion length L v (T) = ^/ctt^/Nq <C 6, where N is the 3d 
density of states. The inequality is satisfied for tempera- 
tures T > T x = To(bk F )- 3 l3 3 / 2 (e F T iPiQ /h) , where e F ,k F 
are Fermi energy and wave number, respectively, and 
l/ r ¥>.o ~ 1-f is the phase relaxation rate at T = Tq = IK 
and /? is a coefficient of order unity. For not too thin sam- 
ples, bk F > 10 , the 3d condition is seen to be met in our 
experiments. 

To gain a first impression of the critical behavior ex- 
hibited by the data we use a fit to the expression 

*(T;Ro)/L = BT» + A. (3) 

Here Lq — e 2 /h is the conductance quantum. The func- 
tional form of the conductivity given by Eq. ((3j> provides 
a good description of all our data for fourteen samples 
in the range 15 kVt < R < 31 fcfi. We fit the data to 
Eq. ([3]) allowing the parameters p, B and A to vary with 
disorder. Figure Q] shows the fit for the sample with a 
value of i?o = 22.67 kVl closest to critical disorder R c 
where A — 0. We note that a positive (negative) value 
of A indicates the dclocalizcd (localized) regime. It fol- 
lows from the scaling relation Eq. (JTJ) that A = at the 
critical point. Fortunately, the sample depicted in Fig. Q] 
with Rq = 22.67 kVl and A = —0.01(1) is close enough 
to criticality so that we can confidently make the iden- 
tification R c = 22.67 fcfi. We note that R c is of order 
unity in units of h/e 2 . Using the same fitting proce- 
dures for the thirteen samples away from criticality (i.e., 
A / 0) we allowed all three parameters, A,B,p to vary 
with Rq. It turns out that B and p do vary with disor- 
der, if only slightly. The scaling property Eq. (JXJ) requires 
the parameters B,p to be independent of disorder. At 
criticality shown in Fig. Q] we determine the parameter 
values (see inset) p = 0.390(3) and B = 0.622(9)if- p . 
We infer from the value of p a critical dynamical expo- 
nent z = 1/p ps 2.5. 

Next we analyse the data in a completely unbiased way 
by using a scaling plot of the scaling function G defined 
in Eq. ([2|). For practical reasons this requires picking a 
value of the critical disorder, R Cl which we take from the 
above analysis as R c = 22.67 k£l. We may now take dif- 
ferent values of the critical exponents v' and v to see how 
well the data points collapse onto single curves on both 
sides of the transition. The best data collapse is found 
for v' = 1.38 on the metallic side and v = 0.77 on the in- 
sulating side. The two panels of Fig. [2]show how well the 
temperature-dependent data for each of the samples in- 
dicated in the legends collapse onto linear scaling curves 
for the two best fit values of v' and v. For convenience, 
we have normalized the axes of each panel to unity using 
the highest temperature value of the conductivity of the 
samples Rq = 21.54 kfl (Rq = 23.77 fcfi) closest to criti- 
cality on the metallic (insulating) side of the transition. 
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FIG. 1: Conductivity as a function of temperature for the 
sample with critical resistance R c — 22.67fcf2 closest to critical 
disorder where A = 0. The solid red line is a fit using Eq. |(3J. 
The parameter A has positive (negative) values for samples 
with Ro < R c (Ro > Rc) 



In the insets to the two panels of Fig. [3J we show 
the dependence of Chi-square (x 2 ) on v 1 in the metallic 
regime and v in the insulating regime. The x 2 calcula- 
tion is performed using statistical weights proportional 
to the inverse of the ordinate values, thereby increasing 
the sensitivity of x 2 to the data sets close to the origin 
and hence further away from criticality. The minima of 
X 2 are clearly defined for p = 0.39 at v 1 = 1.38(11) in 
the metallic regime and v = 0.77(11) in the insulating 
regime. Accordingly, the conductivity exponent s in 3d 
has been experimentally determined to have the value 
s = v' = 1.38(11). 

It is remarkable and somewhat surprising that the 
scaling functions are very well represented by straight 
lines in the regime considered. Thus the conductivity 
on the metallic side can be represented by an expression 
reminiscent of Eq. ([3]) with, however, fixed coefficients 
B,p = l/z, a(T;R )/L = BT 1 '* + a\R - R c \ u ' , and 
a disorder independent coefficient a determined from the 
intercepts of the linear fits in Fig. f2J 

As pointed out in Ref . Q , reasons for the earlier exper- 
iments not agreeing with one another have been traced to 
difficulties in having a system allowing sufficient access 
into the critical region and possessing a well defined crit- 
ical point. In contrast, the current work is done on a sys- 
tem where the critical point can be clearly identified, and 
the critical region is experimentally accessible. Note that 
the number of data sets (14 total) around R c to be kept in 
the scaling analysis of Fig. [2] were determined by compar- 
ing the x 2 fits for different number of data points kept. A 
2 was obtained when the data points were 
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FIG. 2: Collapse of the data for different values of Ro accord- 
ing to the scaling relation given by Eq. ([2]) for v' — 1.38 on the 
metallic side (top panel) and v — 0.77 on the insulating side 
(bottom panel). Here To = 5 K is a reference temperature. 
The insets of each panel show the dependence of \ 2 ° n v> an d 
v for the values of p listed in the legends. The best-fit values 
occur at p = 0.39 with well-defined minima at u 1 — 1.38(11) 
and v = 0.77(11). 



ically estimate the width of the critical region from the 
boundary frequency, w ? = i(£/7)~ z ~ 7(^0 - Rc\/R c ) z , 
where z is equal to u'z on the metallic and vz on the 
insulating side. At T ~ 20K and using 1/r ~ 1Q 3 K, this 
gives |i?o — Rc\/Rc ~ 0.4 on the metallic side and ~ 0.2 
on the insulating side, corresponding to a critical region 
extending from about 15 kQ — 28 fcfi, as shown in Fig. |3] 
This is a sufficiently large experimentally accessible re- 
gion that allows us to obtain the critical exponents quite 
reliably. 

Although the annealing treatment of a sample in the 
critical regime allows one to move the system through 
the transition point in small steps, it is desirable to have 
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FIG. 3: Estimate of the width of the critical regime as a 
function of temperature, using the experimentally obtained 
values of the exponents v' and v. 

in addition a continuous control parameter. We have ob- 
served that an applied magnetic field H shifts the param- 
eter A smoothly. A continuous variation of H therefore 
should allow to access the critical behavior of A in the 
range one order of magnitude smaller than the present 
one. Work in this direction is in progress. 

Our finding of the two distinctly different values for the 
critical exponent v in the metallic and insulating phase, 
respectively, is unexpected. The quality of the data and 
of the fit to the scaling function appears to be so good 
that the difference of v and v' can not be explained by 
uncertainties in the measurement or deviations from the 
scaling form. From the point of view of theory, different 
critical exponents of the correlation length on both sides 
of the transition may indicate a different structure of crit- 
ical modes. We can only speculate that this behavior 
may be related to the system considered here to be ferro- 
magnetic. In ferromagnets we expect spin transport by 
spin wave propagation, which should render it different 
from charge transport. To our knowledge a theoretical 
treatment of this problem is not available at present. 

To summarize, we have studied the metal-insulator 
transition in thin ferromagnetic disordered films in an 
effectively three-dimensional regime experimentally and 
theoretically. We concentrate on the critical regime, for 
which the dynamical scaling behavior at the critical point 
is known to be o{uS) oc a; 1 / 2 . At finite temperatures the 
relevant frequencies uj are given by the temperature T. 
A first fit of the a(T) data to a temperature power law 
plus a T-indcpcndent constant A(Rq) showed an expo- 
nent p ~ 0.4 and serves to determine the critical resis- 
tance R c from the zero of A(Rq). A full scaling analysis 



of the data from a reasonably accessible regime around 
the critical point (±30%) allows to determine the criti- 
cal exponent of the correlation length as v' = 1.38(11) 
on the metallic side and the localization length exponent 
v = 0.77(11) on the insulating side and the dynamical 
critical exponent as z « 2.5. 
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